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ABSTRACT 



Experiments were carried out using an existing tube 
alone, and with the tube connected to a two-dimensional 
wedge nozzle. The range of maximum duration of steady 
reflected pressure from 3.5 to 5 milliseconds was achieved 
through tailored operation for incident shock strengths of 
3.4 and 2.0, using pure Helium and a 70% Helium/30% Nitrogen 
mixture as the driver gas respectfully. 

Spark and continuous light shadowgraph techniques were 
attempted using an optical window at the Mach 4.3 location. 
Results demonstrated that the short duration flow phenomena 
in a shock tunnel can be recorded successfully using 
existing equipment. Calculations showed that the addition 
of a Mach 10 nozzle and 15 m 3 (6' diameter x 15' long) dump 
chamber would provide a useful hypersonic facility for 
instruction and research. 
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I. 



INTRODUCTION 



The simulation of flight conditions using ground based 
facilities has been an essential part of all aircraft 
development. Figure 1 illustrates the various types of 
facilities which have evolved to simulate flight at 
progressively higher speeds and altitudes. Increasing 
energy requirements led to blow-down wind tunnels. The need 
to avoid condensation in nozzle expansions required heaters. 
The very high stagnation enthalpies associated with 
hypersonic flight led to very short duration, intermittent 
techniques. Each of the types shown in Figure 1 was in 
operation more than 20 years ago. A renewed interest in 
hypersonic simulation has been generated by the National 
Aerospace Plane (NASP) program. 

A recent review of facility requirements for hypersonic 
propulsion system testing by Dunn and Lori [Ref. 1] 
concludes that a single ground based facility will not be 
capable of providing the necessary testing to cover the NASP 
flight envelope. The planned use of Computational Fluid 
Dynamic (CFD) codes to bridge the gap left by inadequate 
ground based facilities will make CFD codes an essential 
part of the design process but, before confidence can be 
placed in these codes, fundamental experiments are needed to 
validate the models incorporated for real gas behavior. 
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Ground Based Testing Facilities [Ref. 1] 



Recognizing the need for CFD validation, the feasibility 
of establishing a shock tunnel capability by modifying the 
existing Naval Postgraduate School (NPS) three inch diameter 
shock tube, was examined in the present study. The time 
duration and magnitude of constant pressure available in the 
reflected region of the shock tube were evaluated. The 
concept of tailoring was examined using two types of driver 
gases, pure Helium (He) and a mixture of 70% Helium and 30% 
Nitrogen (N 2 ) , driven into Nitrogen as a test gas. An 
exploratory shock-tunnel nozzle experiment was conducted, 
and the capability of recording short-duration flow 
phenomena using available visualization techniques was 
demonstrated . 

It was concluded that the development of an operational 
hypersonic shock tunnel was feasible and would provide the 
means for follow-on efforts in the area of CFD code 
validation. At the same time, such a facility would provide 
a useful illustration of hypersonic flow to aeronautical and 
space applications engineering students. 

The following material is divided into three sections; 

1) a review of shock tube and shock tunnel ideas, 2) a 
discussion of the experimental investigation including a 
description of the facilities, experimental procedure and 
results; and (3) an evaluation of nozzle throat and dump 
chamber sizing requirements for a hypersonic tunnel. 
Instrumentation modifcations will be required in the 
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development of a hypersonic capability. Appendix A reviews 
the existing facility and provides recommendations for 
upgrading the current hardware and instrumentation. 
Simplified flow models are developed in Appendices B and C 
to enable the selection of nozzle throat diameter and dump 
chamber volume respectively. Experimental pressure and thin 
film traces obtained in the experimental program are 
contained in Appendix D together with predicted results 
obtained using a computer performance code, which is listed 
in Appendix E. 
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II. THEORETICAL CONSIDERATIONS — A REVIEW 



A. BASIC SHOCK TUBE 

A shock tube is an economical short duration test 
facility used in laboratory experiments in gas dynamics, 
aerodynamics and physics. The principles of shock tube 
operation , which are well documented in numerous texts, for 
example, Gaydon and Hurle [Ref. 2] and Owczarek [Ref. 3], 
are illustrated in Figures 2 and 3 [Ref. 4]. As shown in 
Figure 2, the device consists of a long tube divided into a 
high pressure and a low pressure section by a diaphragm. 

When the diaphragm separating the high pressure driver 
section, region 4, and the low pressure section, region 1, 
is ruptured, a series of compression waves coalesce into a 
single shock front (S) which compresses and heats the region 
1 test gas to region 2 conditions. Simultaneously, a series 
of rarefaction (expansion) waves (R) travel into the high 
pressure section, region 4, cooling and expanding the driver 
gas to region 3 conditions. 

Region 3 (the expanded driver gas) and region 2 (the 
compressed test gas) are separated by a contact surface (C) . 
The gases on either side of the contact surface travel with 
the same velocity, and are at the same pressure, but are 
different in temperature. A time history of the process 
(time vs. distance) is shown in Figure 3. 
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Figure 2. Thermodynamic States in the Operation of 
a Shock Tube [Ref. 4] 
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Figure 3. Time History of Events after Diaphragm 
is Broken [Ref. 4] 
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What is shown in Figure 3 can be determined accurately 
using the method of characteristics and shock theory (see 
[Ref. 3]). As the shock front reflects off the driven tube 
endwall, the region 2 gas is effectively shocked a second 
time, creating in region 5 stagnant gas at high pressure and 
temperature. This region may be used for high temperature 
gas dynamic studies or, when coupled to a nozzle, made to 
serve as a reservoir of compressed and heated gas for a 
hypersonic wind tunnel. 

B. PERFORMANCE 

The strength of the generated shock (M s ) is related to 
the pressure ratio (P^/P^ across the diaphragm and the speed 
of sound in the undisturbed gases through the equation given 
by Gaydon [Ref. 2]. 




2viM s — (Yj— 1) a-A 

(Yi+1) (Yj+l) a 4 / 




( 1 ) 



Figures 4 [Ref. 5] and 5 [Ref. 6], which illustrate Equation 
(1), show that the strongest shocks are generated by using 
high pressure ratios and by using driver gases with high 
speeds of sound and low ratios of specific heats. 

As the desired shock strength increases, the required 
P^/P-t increases significantly. To achieve stronger shocks 
several approaches are possible. The most common include 
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Figure 4. Pressure Ratio (P^/P^ Required to 

Generate a Given Shock Strength Using 
Various Gas Mixtures [Ref. 5] 
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Figure 5. Variation of Shock Strength as a Function 
of Speed of Sound 
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Figure 5 (CONTINUED) 
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the use of an area reduction between the driver and driven 
sections [Ref. 5], the use of multiple diaphragms with a 
buffer gas in the driver [Ref. 5], and heating the driver 
gas through combustion or electrical discharge techniques 
[Ref. 6]. Figures 6 and 7 [Ref. 7] show the relationship 
between shock strength (M s ) and reflected region 5 
conditions for equilibrium air. 

C. SHOCK TUNNELS 

The reflected region 5 pressure and temperature may be 
used to supply a supersonic or hypersonic nozzle if the 
nozzle throat is made sufficiently small to ensure full 
reflection of the incident shock. The evolution of this 
arrangement is shown schematically in Figure 8. To maximize 
test time and minimize the nozzle starting time, often a 
diaphragm separates the nozzle and dump chamber from the 
driven end of the shock tube. The dump chamber is sized and 
evacuated to a level which ensures that, when the tube is 
fired, the pressure in the chamber remains below the static 
pressure at the nozzle exit for the duration of steady flow. 
[Ref 8 : pp . 624-640] 

For hypersonic applications a conical nozzle is often 
used. Table 1 provides calculated test conditions for 
various flow expansions using ideal gas theory. The 
computer program used for the calculations is described in 
Appendix E. Since, in an isentropic expansion, static 
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Figure 6. Reflected Pressure as a Function of Shock 
Strength for Equilibrium Air [Ref. 7] 
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Figure 8. Evolution of a Shock Tunnel 
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TABLE 1 



THEORETICAL PREDICTION OF TEST CONDITIONS FOR A HYPERSONIC 
NOZZLE (M e = 10) DRIVEN FROM THE REFLECTED REGION 

OF A SHOCK TUBE 



Driver 




P 5 


t 5 


P e 


T e 


^■start 


Gas 


M s 


(Atm) 


(°K) 


( mm Hg ) 


(°K) 


(msec) 




” 2.4 


58 


973 


1.02 


50 


. 041 




3.7 


32 


2000 


.55 


94 


. 029 


Helium 


4.2 


23 


2460 


.41 


116 


. 026 




4 . 5 


19 


2700 


.33 


131 


. 024 




~ 1.8 


23 


632 


.4 


30 


. 051 


Nitrogen 


2.4 


9.7 


971 


. 169 


50 


. 041 




2 . 6 


6.7 


1110 


. 118 


52 


. 038 




2 . 8 


5.3 


1200 


.09 


56 


. 037 



temperature decreases as the Mach number (of the test 
section) increases, liquefaction of the test gas becomes an 
issue. Since nitrogen closely matches the aerodynamic 
properties of air, but liquefies at a lower temperature, it 
is preferred for hypersonic facilities in which simulated 
real gas effects are not required. This is illustrated in 
Figure 9 which also shows the interesting results of wind 
tunnel experience [Ref. 9]. In practice, condensation is 
observed to occur at temperatures which are lower than the 
saturation temperature, thereby reducing somewhat the 
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Figure 9. Experimentally Determined Liquif ication 
Properties of Air and Nitrogen [Ref. 9] 
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stagnation temperature requirements to reach a given Mach 
number. 

Given the extremely short test time available in a shock 
tube, any coupling to a hypersonic nozzle must be done in 
such a way so as to minimize nozzle starting transients and 
provide the maximum available steady flow in the test sec- 
tion. The starting time, defined as the time to approach 
steady flow through a hypersonic conical nozzle ( 7 X = 7/5) , 
is given by Smith [Ref. 8] as 



. 76Rj_ 

^start 8 a c 



-IM 



1/2 



(1 



e 



TF (M e 



3/2 



- 1 ) + 



2? (M e 



5/2 



-D] 



(2) 



A reduction in starting time may be achieved by using a 
large nozzle angle (ft = nozzle half angle) , a small throat 
radius (R t ) and a high stagnation temperature in the 
reflected region. [Ref. 8:pp. 625-640] The variation in 
starting time with increasing incident shock strength for a 
conical nozzle with ft = 30 degrees is shown in Table 1. 

After the unsteady starting process, the useful test 
flow will continue as long as a sufficiently high pressure 
ratio is maintained across the nozzle. The steady flow 
process will be a function of only the upstream reflected 
stagnation conditions as long as the dump chamber has 
sufficient volume. Ideally, uniform reservoir stagnation 
conditions will end when the reflected region conditions are 
changed by a wave emanating from the interaction of the 
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contact surface and the head or tail of the reflected expan- 
sion wave (Figure 2). The useful test time in the nozzle 
ends when the reflected pressure drops below the level 
required to maintain the static pressure at the nozzle exit 
above the back pressure in the dump chamber. Thus the 
useful test time will depend on the nozzle throat size, 
which governs the flow rate out of the tube, and the size of 
the dump chamber, which determines the rate at which the 
back pressure builds up. 

D. TAILORED OPERATION 
1 . Tailoring 

The fundamental limitation of a shock tube is the 
available test time. Considerable effort has gone into 
trying to extend shock tube test times [Ref. 10]. Test 
times can be lengthened by as much as a factor of ten by 
"tailoring, " which means adjusting conditions so that the 
ratio of the speed of sound and specific heat behind the 
reflected shock matches those in the region behind the 
contact surface. When this tailored condition is present, 
the reflected shock passes through the contact surface 
bringing it to rest. A Mach wave is reflected forwards off 
the contact surface and generates no additional disturbances 
as it passes through region 5. In the untailored condition, 
the wave reflected forward from the contact surface is not a 
Mach wave but a rarefaction wave when aj < a 3 , and a shock 
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wave when a : > a 3 . The reflected waves cause conditions in 
region 5 to change. Experimental evidence of the effect of 
the wave interactions obtained by Pennelegion [Ref. 10] are 
shown in Figure 10. Pennelegion ' s research fully explored 
tailoring analytically and provided supporting experimental 
evidence. His findings and experimental techniques were 
followed in the experiments undertaken in the present study. 
The tailored interface technique requires that the following 
equation be satisfied [Ref.l0:p. 4]; 



2i 

y 4 



Y-j+1 P5 

1+ 

V 1 P 5 

- 1 + ( iV } ( p^ _1) - 

y 4 2 



1/2 



( 3 ) 



A unique tailored condition (shock Mach number) , exists for 
every combination of driver and driven gas. The duration of 
constant pressure in region 5 lasts until the rarefaction 
wave reflects from the driver endplate and returns to meet 
the contact surface. For low incident shock strengths it is 
the head of the rarefaction wave, and for high incident 
shock strengths it is the tail of the rarefaction wave which 
terminates the constant pressure. 

The tailored interface technique for producing 
extended periods of constant reflected pressure and 
temperature was demonstrated experimentally for Mach numbers 
below 3.4. When the technique was applied to shocks of Mach 
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Figure 10. Experimental Observation of the Three Types 
of Wave Interactions in the Reflected Region 
of a Shock tube [Ref. 10] 
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number greater than 3.8, there was a significant drop in 
test time due to boundary layer interaction. In the under- 
tailored condition the reflected conditions showed an 
initial fall in temperature and pressure, followed by a slow 
rise in pressure and continued drop in temperature. The 
overtailored condition revealed a pressure rise from the 
reflection from the contact surface followed by a rapid drop 
in temperature. In all cases with strong incident shock 
waves, the real gas effects caused significant variation in 
the reflected region conditions. [Ref. 10] 

2 . Mixing of Driver Gases (He and N-, ) 

The application of tailoring to increase test time 
in a hypersonic shock tunnel allows more accurate 
measurements of force and pressure to be made on test 
models. As desired test section Mach number increases there 
is a significant decrease in working section static 
pressure. It is therefore desirable to have reflected 
region pressures as steady and as high as possible. Pressure 
recovery, defined as the ratio of the initial step in the 
reflected pressure (P 5 ) to driven pressure (PJ , may be 
increased by increasing the speed of sound of the driver if 
some test time is sacrificed. The combination of a driver 
of Helium (He) and Nitrogen (N 2 ) can provide an improvement 
in test time and pressure recovery over a driver of N 2 
alone. Experimental evidence of this is shown in Figures 11 
and 12 from Pennelegion [Ref. 10] . From the same study, 
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Figure 12. Experimental Pressure Recovery for Different 
Driver Gases [Ref. 10] 



23 



Figure 13 and Table 2 show tailoring using various mixtures 
of He and N 2 to generate shocks in the range 1 < M s < 3.4. 

The tailoring shock Mach number may be selected to provide 
the required stagnation temperature to prevent liquefaction 
and maximize test time for hypersonic nozzles with various 
Mach numbers. Pennelgion [Ref. 10], provided the method for 
calculating tailored conditions for He, N 2 mixtures. 




0-3 0 6 0 * 7 0 - 8 0-9 1 0 

helium fraction 



Figure 13. Variation of Tailored Mach Number with 
Helium Fraction [Ref. 10] 
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TABLE 2 



EXPERIMENTALLY TAILORED MACH NUMBERS FOR DRIVER 
GAS MIXTURES [REF. 10] 



He 

Fraction 


Mte 


Pa/ Pi 


Time Steady 
Pressure (msec) 


t 5 

(°K) 


0.7 


2.0 


24 


6.5 


714 


0.8 


2.7 


60 


6.0 


1096 


0.9 


2.8 


70 


4.7 


1159 


1.0 


3.4 


124 


4.3 


1553 
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III. SHOCK TUBE 



In this section the existing shock tube and instrumenta- 
tion, and the modifications made during the course of the 
present experiments, are described. 

A. STRUCTURE 

The three inch diameter shock tube at the Naval 
Postgraduate School Aeropropulsion Laboratories is shown in 
Figure 14. A detailed description is given in Penaranda 
[Ref. 4]. The tube, constructed of highly polished 303 
stainless steel consists of a driver section 7.75 feet long, 
shown in detail in Figure 15, and a driven section 18.9 feet 
long . 

B. DIAPHRAGM CHAMBER 

The tube incorporates a double diaphragm arrangement to 
maintain precise and repeatable diaphragm bursting pressure 
ratios (P^/PJ . The tube currently is fitted with diaphragms 
designed to burst using a constant driver pressure of 600 
psi. Scored brass diaphragms designed to rupture at 400 psi 
are inserted. Variation in shock speed is obtained through 
adjustment of the downstream pressure (P x ) . To obtain the 
desired pressure ratio, the 2.5 inch long diaphragm chamber 
(Figures 15 and 16) and the driver chamber are simul- 
taneously filled to one half of the desired P4 pressure (300 
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Figure 15. Shock Tube Driver Section [Ref. 4] 
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Figure 16. Dual Diaphragm Chamber (Unbolted) 

psi) . A diverting valve is closed to isolate the diaphragm 
chamber at 300 psi. The driver section pressure is then 
raised to the desired P 4 pressure (600 psi) . To fire the 
tube, a diverting valve is opened and some of the driver gas 
at 600 psi is dumped into the diaphragm chamber causing the 
first diaphragm to break. The second diaphragm breaks im- 
mediately, as the pressure on the driven face of the 
diaphragm rapidly falls. 
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The dual diaphragm arrangement minimizes bursting 
pressure fluctuations due to inconsistencies in the material 
of the diaphragms. By applying a pressure of 200 psi above 
the design bursting pressure of the first diaphragm, rupture 
of the two diaphragms is ensured with a precise value of P A . 
In shock tubes using a single diaphragm arrangement, precise 
bursting pressures cannot be repeated unless the complexity 
of a diaphragm puncturing device is used. Without a 
puncture device, the bursting pressure fluctuations remain a 
function of the uniformity of the diaphragm material. 

C. NOZZLE ASSEMBLY 

An existing two dimensional stainless steel nozzle 
designed for an exit Mach number of 5.72 (A/A* = 43.7) and a 
rectangular dump chamber to which it bolted were available 
for the present study. The assembly is shown in Figure 17. 
For the present experiments two 1" diameter window ports 
were installed to allow optical measurement techniques to be 
exercised, and a model support was constructed to position a 
3/8" axisymmetric model on the nozzle centerline in line 
with the window, where the Mach number was about 4.3. The 
final assembly is shown in Figure 18. The dimensions of the 
nozzle are given in Figure 19 and the design of the window 
installation is shown in Figure 20. Views of the nozzle and 
model support are shown in Figure 21. Since the nozzle and 
dump tank existed (from a previous gas dynamic laser 
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Figure 17. Two Dimensional Nozzle and Dump Chamber 







Figure 18. Assembled Dump Chamber, Model Support 
and Nozzle 




Figure 19. Dimensional Drawing of the Nozzle 
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Figure 20. Optical Window Installation 




(a) Model, Model Support and Nozzle 




(b) Nozzle Exit 

Figure 21. Views of the Nozzle and Model Support 
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investigation) , location and sizing of the window as well as 
the maximum flow duration were limited by the existing 
hardware . 

D. INSTRUMENTATION 
1 . Shock Speed 

Shock speed was determined by measuring the time for 
the shock to travel between thin film gauges located at 
positions A and B and from B to the pressure transducer at 
C (Figure 14) . The locally manufactured thin film gauges 
shown in Figure 22, consist of a 0.25" diameter pyrex rod 
coated with a 1/32" wide thin film of platinum giving a cold 
resistance of 100 ohms. [Ref. 4:p. 23] As the region of 
hot gas immediately behind the shock front passes over the 
platinum, its temperature and resistance change and a 
voltage pulse is generated. In the present work the signals 
from the thin films were amplified and recorded for 
analysis. The thin films at location A also served as a 
trigger for the digital counter and shadowgraph spark 
source . 

Shock speed was determined in different ways. The 
pre-existing techniques involved using oscilloscope records 
of thin film and pressure transducer outputs, and also the 
output of a digital counter triggered on and off by the thin 
films at A to B respectively. In the present work, in runs 
1 through 49 shock speeds were determined from a Polaroid 
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Figure 22. Thin Film Instrumentation Assembly 
[Ref. 4] 
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photograph of the oscilloscope, such as is shown in Figure 
23, and from the digital counter. Both techniques were 
identical to those used in the late 1960's. 

In runs 50 -85, the outputs of the three traces were 
recorded digitally on three of the available four channels 
of a Scientific Alanta SD380 Spectrum Analyzer. Data were 
recorded at .019 and .039 msec intervals depending upon the 
time scale selected. An example is shown in Figure 24. The 
shock speeds were determined during postprocessing of the 
data . 

2 . Driver Tube Pressure 

A dial gauge with 10 psi pressure increments was 
used to determine when the constant 600 psi driver pressure 
was reached. A dead weight tester was used to calibrate the 
gauge . 

3 . Driven Tube Pressure 

Two calibrated Wallace and Tiernan Bourdon tube 
gauges were used to determine the driven tube pressures. In 
the range from 0 to 50mm Hg, a gauge with increments of .2 
mm Hg was used. For higher driven pressures, a second gauge 
with a range from 0 to 35 psi in increments of .05 psi was 
used . 

4 . Driven Tube Temperature 

A Type K (chromel-alumel) thermocouple, connected to 
a Daigger Scientific Digital Panel Meter using matched 
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Figure 23. Typical Polaroid Photograph of the 
Thin Film Response from Location B 
and the Output of the Original 
Pressure Transducer 
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Figure 24. A Typical Recording Using the SD380 Signal Analyzer 



extension wires, was mounted in the driven end of the tube. 
The digital readout was recorded prior to every run. 

5 . Reflected Pressure and Test Duration 

A piezoelectric pressure transducer, Endevco model 
2501, a charge amplifier, Endevco model 2700, and a 
Tektronics 551 oscilloscope were initially used to record 
the transient region 5 pressure. The lack of resolution and 
strange qualitative behavior observed in the results (Figure 
23), led to the mounting of an improved transducer. 

A PCB Piezontonics , Inc. Model 124A pressure trans- 
ducer connected to a model 504A dual gain charge amplifier 
were used in runs 50 through 75. The system was calibrated 
using a standard static calibration feature described in the 
manufacturer's manuals [Refs. 11,12]. The transducer 
provided a response of 240 psi/volts. Experimental traces 
obtained are given in Appendix D. 

6 . Oscilloscopes 

Two Tektronics oscilloscopes, models 549 and 551, 
fitted with Polaroid cameras were used for runs 1 through 
49. The capability of the scopes are fully described in the 
associated manuals [Refs. 13,14]. Time measurements to 
determine shock speed, and voltage levels to determine 
pressure, were determined by overlaying a scale on the 
Polaroid photographs. 
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7 . Signal Analyzer 

A Scientific Alanta SD380 four channel signal 
analyzer (Figure 25) equipped with a hard disk and plotter 
was used to record and process data from runs 50 through 85. 
The independent channel trigger delays (±) , adjustable 
scales and post-processing capabilities of the digitized 
data (in contrast to photographing oscilloscopes), is 
evidence of the advances in measurement capabilities over 
the past 20 years. A full description of capabilities of 
the analyzer is given in the instrument manual. [Ref. 15] 

E. OPTICAL MEASUREMENT SYSTEMS 

A simple shadowgraph arrangement was generated from the 
components of a 6" schlieren system. The optical layout is 
shown in Figure 26 [Ref. 16]. The spark source was used 
with Polaroid film in a camera back attached to the nozzle 
window. The continuous source was used with a Hycam movie 
camera . 

1 . Light Sources 

a. High Pressure Mercury Arc Lamp 

A high pressure mercury arc lamp (Per Type 112- 
5145) constant light source was used with a six inch 
achromatic convex lens with a focal length of 60 inches. The 
arrangement, mounted on an adjustable bench (Figure 27), 
provided a parallel light beam through the nozzle windows. 
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Figure 25. Scientific Alanta SD380 Four 
Channel Signal Analyzer with 
Hard Disk Installed 
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Figure 26 . Optical Layout Sketch 
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Figure 28. Spark Source with High Voltage 
Power Supply 
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The beam was adjusted by the procedures described in the 
associated manual [Ref. 16]. 

b. NASA 229-D Spark Light Source 

A NASA 229-D spark light source (Figure 28) 
provided a nanosecond spark for shadowgraph pictures taken 
with Polaroid film. The spark was triggered from a thin 
film connected through the delay and voltage circuits shown 
in Figure 29. The timing for the spark source was 
determined experimentally. 

2 . Camera Systems 

a. Hycam Movie Camera 

A Hycam movie camera [Ref. 17] with a framing 
rate on the order of 6000 pictures per second at runout of a 
400 ft roll of film, was used. The camera was equipped with 
a Nikon 55mm lens and ASA 250 black and white film. A time 
reference was provided on the edge of the film by a strobe 
light pulsing at 1000 pulses per second. The camera was 
triggered manually. The camera was mounted on a second 
adjustable bench as shown in Figure 30. 

b. Spark Shadowgraphs 

Spark shadowgraphs were taken using a Polaroid 
Land Film Holder with Type 52 ASA 400 Black and White film, 
held directly behind the optical windows. To eliminate 
stray light, the mounting arrangement was contained in a 
cardboard box, and tests were conducted after dark. 
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Figure 29. Spark Timing Delay Circuit and High 
Voltage Pulse Generator 



Figure 30. Hycam Movie Camera with Separate 
Pulse Generator 
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IV. EXPERIMENTAL PROCEDURE 



A. PREPARATION 

Prior to testing, with the driven endplate removed and 
the diaphragm chamber opened, a cloth lightly dampened with 
Acetone was pulled through the driver and driven sections. 
O-ring grooves of the opened sections where thoroughly 
cleaned and a thin coating of vacuum grease was applied to 
the rubber O-rings. Diaphragms were wiped clean of visible 
dirt and the system reassembled. The cleaning procedure was 
repeated after approximately 15 runs. Using a roughing 
pump, a clean-tube vacuum pressure of 3 mm Hg was routinely 
achieved. 

B. LOADING THE DRIVEN GAS 

Prior to loading to the desired test pressure (P a ) , the 
driven section was purged. To purge, the tube was evacuated 
to 3 mm Hg and then isolated from the roughing pump. 

Nitrogen was introduced to a pressure 1 to 4 psi. The tube 
was then evacuated a second time and repressurized to the 
desired test condition. 

C. LOADING THE DRIVER GAS AND FIRING 

The driver assembly is designed to operate at 600 psi 
with single or dual gas combinations. To purge the driver 
section and diaphragm chamber, the system was pressurized to 
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approximately 30 psi and then vented to atmospheric pres- 
sure. The system was then pressurized to 600 psi with the 
diverter valve to the diaphragm chamber being closed at 300 
psi. When the 70% He and 30% N z driver mixture was used, 
the system was pressurized to 420 psi with He and then to 
600 psi with N 2 . 

The tube was fired by opening the diverter valve to the 
diaphragm chamber. Some seepage from the driver section was 
evident. To minimize this effect the tube was fired 
immediately after reaching the 600 psi level. Replacement 
of 0-rings would minimize this problem. Seepage into the 
driven section was not evident. 
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V. EXPERIMENTAL RESULTS AND DISCUSSION 



A series of experiments was conducted to gain experience 
with the shock tube. The study consisted of approximately 
85 shots involving the use of He and a 70% He/30% N 2 mixture 
as driver gases, driven into a N 2 test gas. Studies with 
the shock tube assembled with the pressure transducer 
mounted in the end of the driven tube, were concentrated on 
investigating the concept of tailoring to maximize the 
magnitude and duration of steady pressure in the reflected 
region. While more than 75 shots were made in this 
configuration, the replacement of the pressure transducer 
and postprocessing of data taken with the SD380 greatly 
enhanced the results. For this reason, only data taken after 
run 50 using the SD380 are presented. Studies with the 
nozzle attached at the end of the driven tube were 
exploratory attempts to obtain flow visualization in a shock 
tunnel configuration. The results of the two series of 
experiments are discussed in the following sections. 

A. TAILORING 

Using the nomenclature of Glass [Ref. 18], a tailored 
condition is defined as the condition when a Mach wave is 
reflected off the contact surface and causes no change in 
the reflected pressure (P 5 ) . This can be identified 



49 



experimentally (Figure 10) as the condition when the ratio 
of the initial reflected pressure (P 5 ) and the pressure 
after the wave interaction at the contact surface reflects 
off the endplate (P 8 ) being equal to one. Figure 31 shows 
examples of results in the present study for the end wall 
pressure untailored and tailored conditions using He with an 
incident shock strength, M s = 3.4. Similar traces for near 
to a tailored mixture (.7 He and .3 N 2 ) with M s = 2.0 are 
shown in Figure 32. These results are in complete agreement 
with Pennelegion [Ref. 10]. 

B. TIME OF STEADY PRESSURE 

It was experimentally determined in the present study 
that a tailored mixture as a driver gas produced a longer 
duration of steady pressure than that of a pure He driver. 
The time of steady pressure for the existing NPS shock tube 
using He as a driver was approximately 3.6 milliseconds 
while the mixture extended the test time duration to 5 
milliseconds. A comparisons of the test times obtained in 
the present study with those of Pennelegion [Ref. 10] is 
shown in Figure 33. The trends are seen to be similar. 

Since the work of Pennelegion [Ref. 10] was performed at a 
constant driver pressure of 1000 psi in a N.P.L. 2" shock 
tube with a 10 ft driver and a 12 foot driven section some 
differences in the magnitudes of the test times are to be 
expected. 
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RUN 58 (MACH = 3.8) 
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Figure 31 




(b) Tailored 



Experimentally Determined Endwall Pressure 
Traces Using Helium as a Driver Gas 
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Figure 32. Experimentally Determined Endwall Pressure 
Traces Using a Mixture of 70% Helium and 
30% Nitrogen as a Driver Gas 
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TIME OF STEADY PRESSURE (MSEC) 




MACH ft 



Figure 33. Comparison of Steady Pressure Time Achieved 
in the NPS 3 Inch Shock Tube with Experi- 
mental Results of Pennelegion [Ref. 10] 
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For low strength shocks the head of the expansion wave 
causes an early degradation of the reflected pressure in the 
NPS tube. The longer driver section in the N.P.L. tube is 
responsible for the difference. 

C. PRESSURE RECOVERY 

Figure 34 shows the fraction of initial driver pressure 
available in the reflected region as a function of shock 
Mach number. The data show that as shock strength increases 
pressure recovery decreases for both He and mixed driver 
gases. The data are in good agreement with the results of 
Pennelegion [Ref. 10]. The measured values of the reflected 
shock pressure, which are the pressures available to drive a 
hypersonic tunnel, are given in Tables 3 and 4. 

D. REFLECTED TEMPERATURE 

The reflected temperature was not determined 
experimentally. Instead, shown in Table 2 are values of 
stagnation temperatures calculated for shock strengths in 
the range from 2.0 to 3.4. These values were taken from the 
tabulated solutions (based on equilibrium chemistry) of 
Bernstein [Ref. 19:pp. 32-43], 

E. SHOCK ATTENUATION 

A representative illustration of measured shock 
attenuation for both a pure gas driver and the driver gas 
mixture is shown in Figure 35. Velocity can be inferred 
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Figure 34. Comparison of Reflected Pressure Measurements 
in the NPS 3 Inch Shock Tube with Experimen- 
tal Results of Pennelegion [Ref. 10] 
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TABLE 3 



EXPERIMENTALLY DETERMINED REFLECTED PRESSURE MAGNITUDE 
AND DURATION AND CALCULATED TEMPERATURE (HE DRIVER/N 2 DRIVEN) 



Run 






P 5 




Time of Ste 


# 


M s 


P*/Pi 


(psi) 


P 8 /P 5 


P 5 (msec) 


50 


3 . 08 


994 


150 


1.79 


1.3 


51 


2.92 


353 


235 


1.23 


1.2 


52 


4 . 39 


414 


207 


1.38 


1.9 


53 


3 .32 


93 


312 


1.42 


1.5 


57 


4 . 52 


500 


183 


1.25 


0.3 


58 


3 .79 


207 


304 


1.23 


0 . 3 


59 


3 . 66 


148 


386 


1.13 


1.6 


60 


3 . 36 


169 


264 


1.13 


3 . 6 


61 


3 .72 


167 


317 


1.23 


0.5 


64 


3.42 


164 






— 


65 


3 . 65 


158 






— 


66 




184 








67 


3.77 


173 








68 


3.81 


400 








74 


4 . 50 


667 








75 


4.10* 


1240 









* Estimated Value 



TABLE 4 



EXPERIMENTALLY 


DETERMINED 


PRESSURE 


MAGNITUDE AND DURATION 


AND 


CALCULATED 


TEMPERATURE 


(70% HE, 


30% N 2 


DRIVER/N 2 DRIVEN) 


Run 






P 5 




Time of Steady 


# 


M s 


P</Pi 


(psi) 


P 8 /P 5 


P 5 (msec) 


62 


2.40 


50 


317 


1.10 


2 . 8 


63 


2 . 66 


51 


310 


1.23 


1.5 


69 


2.47 


50 


310 


1.09 


2 . 5 


70 


2.25 


50 


305 


1.11 


3.3 


71 


2 . 03 


24 


427 


1.02 


4.9 


72 


2 . 03 


22 


422 


1.08 


5 


73 


1.99 


21 


410 


1.13 


2 


77 












78 




155 








79 


3.0* 


171 








80 


3.4* 


638 









* Estimated Value 
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from the slope of the lines between the thin films A and B 
and between B and the pressure transducer at C. An 
increasing slope indicates a slower speed. The general 
trend is that for a given gas mixture, the stronger the 
shock strength the smaller the attenuation between B and C. 
Reflected velocity is of the order of one half the incident 
velocity. Tabulated results for all runs are found in 
Appendix D. 

F. OPTICAL MEASUREMENTS 

Before attempting to obtain photographs during the 
firing of the shock tube, a small freejet from a 1/4" 
flattened copper tube, powered by shop air, was installed in 
the model location, the spark light was adjusted using a 
neutral filter and a successful shadowgraph was taken. 

1 . Shadowgraph Using Hvcam Camera 

The Hycam Movie camera was triggered manually either 
one or two seconds before the diaphragm was ruptured. The 
400 ft rolls of ASA 250 speed film provided approximately 
one second for the camera to accelerate, and a three second 
window to capture the event. It was found that all the 
available light was needed to obtain a correct exposure and 
the camera was focused on the outside surface of the window. 
The shadowgraph was made visible by the faint film which by 
then had deposited on the inside of the window. The result 
was successful. 
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With the blunt body shown in Figure 21(a) installed, 
the nozzle starting Mach wave (Figure 36(a) and (b) ) , and 
the establishment of steady flow with an oblique shock 
visible (Figure 36(c)), were successfully recorded. Steady 
flow (started) conditions lasted a brief time until dump 
chamber back pressure unstarted the nozzle. Incident shock 
strength was estimated to be M s = 3.8 providing a peak 
pressure in the reflected region of approximately 460 psi. 
The actual shock speed was not measured during the test 
since the starting of the camera caused premature triggering 
of the spectrum analyzer. However the initial conditions, 
and therefore the reflected pressure, were the same as in 
run 58. A complete record of this run given in Appendix D. 

2 . Spark Shadowgraph 

Before firing the shock tube, the correct exposure 
for the expected shadowgraph was verified by opening the 
film and firing the spark source. The result is shown in 
Figure 37. It was subsequently found that when the shock 
tube was fired, the thin film signal failed to trigger the 
spark, and thus any exposure on the film was the result of 
light produced from within the tube. The results of 
exposing the film without discharging the spark as the shock 
Mach number was progressively increased, are shown in 
sections (a) to (e) of Figure 38. A faint glow appeared 
first at a M s = 3.6. As M s was increased for both a pure He 
driver and mixed driver gases, the luminosity increased. 
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(a) Nozzle Starting Shock Wave 




(b) Nozzle Starting Shock Wave 

Oblique Shock Wave 







(c) Oblique Shock 

Figure 36. Results Obtained Using Hycam Camera 
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Figure 37. Spark Shadowgraph with No Flow Present 
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Film Exposure without a Spark Source 
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Figure 38. (CONTINUED) 
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Figure 38. (CONTINUED) 




The source of the light requires further investigation, but 
it is clear that some contamination must be present since 
the driver and driven gases, if pure, would not be excited at 
the temperature produced in the firing process. 



65 



VI. SIZING CONSIDERATIONS FOR A HYPERSONIC TUNNEL 



The experimental program showed that the basic shock 
tube gave reliable and repeatable results and, following the 
upgrading recommendations detailed in Appendix A, it could 
well serve as a driver for a hypersonic tunnel. Having 
experimentally defined the stagnation temperature and 
pressure available to drive such a facility, it only 
remained to arrive at a nozzle and a dump chamber 
configuration suitably matched to the tube. Simplified 
models were therefore developed to, first, calculate the 
effect of nozzle throat diameter on the available steady 
flow time and, second, to determine a suitable dump chamber 
size. The two calculations are detailed in Appendix B and 
Appendix C respectively. 

A. NOZZLE THROAT SIZE 

The nozzle throat diameter must be sized such that the 
time for which conditions in the reflected region are steady 
must be greater than the time required for the compressed 
mass of gas to exhaust through the nozzle. The model 
described in Appendix B estimates the time in relation to 
the falling pressure in the reflected region to be given by 




(B. 10) 



P 



t 



L 5 



P, 
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If it is assumed the supply pressure can drop to half its 
initial value, the time required can be calculated using 
Equation (B.10) for different driver gases and shock 
strengths. Through the ratio m 0 /m 0 , the calculated time is 
implicitly a function of the throat diameter. Results of 
such calculations are shown in Figure 39(a) and (b) . The 
four curves shown illustrate the dependence of useful time 
on the shock strength and throat diameter. Clearly, the 
shock strength is less important than the throat area. 

Based on these results, a throat diameter 6.5mm (.025 
inches) was selected. This choice made best use of the 
steady flow time which was experimentally determined to be 
from 3.5 to 5 milliseconds using tailored conditions. For 
the 6.5mm throat, Equation (B.10) was used to calculate the 
time for the reflected shock region pressure to decay to 
various fractions of its initial magnitude. The results are 
shown in Figure 40. Clearly, for even the strongest shock, 
the time for the pressure to decay to 80% of its initial 
value is longer than its experimentally determined tailored 
test time of 5 msecs. 

B. DUMP CHAMBER SIZE 

The dump chamber requires a capacity such that the back 
pressure remains below the nozzle exit pressure for the 
duration of steady conditions in the driven tube. The 
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(a) Helium into Nitrogen 



Figure 39. Time for Reflected Region Pressure to Drop 
to Half its Initial Value 
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TIME TO RERCH ONE-nnLF P5 (SEC) 




(b) Nitrogen into Nitrogen 



Figure 39 (CONTINUED) 
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Calculated Time of Decay of Reflected 
Region Pressure for a 6.5mm Nozzle Throat 
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analysis in Appenix C gave the required ratio of dump 



chamber volume to reflected shock reservoir volume as 





V 



D 



V 



0 



Examination reveals that the parameters which most influence 
the required dump chamber volume are the incident shock 
strength (which determines the reflected reservoir 
conditions) , the desired test section Mach number, the 
driver gas selected, and the specified reflected region 
pressure (P t /P t o) at which the flow breaks down. Since the 
pressure recovery is a function of which driver gas is 
selected, calculations were made for both nitrogen and 
helium. The results are shown in Figure 41. A mixture of 
the driver gases would be expected to fall between these two 
curves. Consistent with the finding above, the criterion 
that the dump chamber pressure becomes equal to the nozzle 
exit pressure when 20% of the stagnation pressure was 
exhausted, was taken to determine the necessary dump chamber 
volume. Shown in the figure are results obtained assuming 
the unknown temperature in the dump chamber to have its 
extreme maximum and minimum values; namely, initial 
reflected shock temperature and initial dump chamber 
temperature. From an examination of the results, a 
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Figure 41. Dump Chamber Volume Requirements 
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conservative selection would be 15 cubic meters. This is a 
vessel approximately 6 feet in diameter and 15 feet long. 

C. HIGH ENTHALPY MODIFICATION 

Higher enthalpy flows can be generated using the free- 
piston arrangement shown in Figure 1. Stalker [Refs. 20,21] 
has carried out extensive experimental work in this area. A 
preliminary analysis shows that the NPS shock tube has the 
same internal diameter as Stalker's 1st generation free- 
piston tunnel, but has a thinner wall by 1/2". Structural 
limitations would prevent the full duplication of Stalker's 
tunnel but, with the published information on the design of 
free piston drivers, it should be possible to extend the 
present tube to an increased enthalpy level at which real 
gas effects of interest can be studied. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



The present study showed that: 

The existing shock-tube is straightforward and 
relatively inexpensive to operate, and gives 
repeatable test conditions; 

The existing shock tube can be used to drive a 
hypersonic shock tunnel to a Mach number of 12 ; 

Tailored operation can provide uniform stagnation 
conditions for periods of approximately 3.4 msec 
using Helium, and approximately 5 msec using a .7 
He/.3N 2 mixture, as the driver gas; 

It is feasible to obtain both high speed movies and 
short duration spark flow visualization records of 
the test section flow field using currently available 
equipment ; 

Contamination in the tube can lead to unwanted 
radiation (which must clearly be eliminated if real 
gas effects are the purpose of the experiment) . 

It is recommended that: 

The existing tube be serviced and its operating 
instrumentation be upgraded, as outlined in Appendix A 

A Mach 10 nozzle, 15 cubic meters dump chamber and 
model support be designed and procured for operation 
with the current tube; 

- Consideration be given to the addition of a free- 
piston tunnel to provide high enthalpy simulation. 

Operation of a shock tunnel with the existing tube would 

provide a convenient and inexpensive facility for use in 

laboratory instruction. It would also provide a tool for 

thesis studies aimed at hypersonic CFD code validation. 
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Extension to a free-piston driver would provide an advanced 
tool for hypersonic research. 
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APPENDIX A 



APPRAISAL OF THE PRESENT SHOCK TUBE 



The NPS shock tube has been in operation since the late 
1960's. This appendix provides a description of recommended 
improvements/enhancements of the existing system, which 
result from the experience gained in the reported program of 
tests . 

A. O-RINGS 

The planned relocation of the facility to the new Gas 
Dynamics Laboratory requires complete disassembly of the 
shock tube. Many O-rings are old and worn and they should 
be replaced at that time. Table 5 details the O-ring 
requirements . 

B. THIN FILM GAUGES 

The unavailability of replacement thin film gauges used 
to measure shock speed is unacceptable. Only the two in- 
stalled gauges remain operational. The options are either 
to manufacture thin films in-house as detailed in Penaranda 
[Ref. 4] or to procure them commercially through the supply 
system. 
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TABLE 5 



O-RING REPLACEMENT GUIDE 





Number 


Width of 


Depth of 


Inside 


Outside 


Location 


Required 


Groove 


Groove 


Diameter 


Diameter 






(inches) 


(inches) 


(inches) 


(inches) 


Driver Tube 


1 


.100 


.157 




3. 502+ . 005 






.106 


.16 3 






Spaced 


2 


.100 


.157 




4.000 


Diaphragm 




.106 


.16 3 




4.005 




2 


.100 


.157 




3.500 






.106 


.16 3 




3.505 


Valve Body 


1 


.101 


.157 


3.2485 








.105 


.163 


3.2495 




Tubing 


3 


.101 


.157 


3.2485 




Sections 




.105 


.105 


3.2495 




Instrumentation 5 


.101 


.157 


3.2485 




Section 




.105 


'.16 3 


3.2495 




Driver End Cap 1 


.101 


.157 


3.2470 








.105 


.163 


3.2460 
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C. PRESSURE TRANSDUCERS 

1 . Piezoelectric Transducer 

The two Endevco Model 2501 transducers with Model 
2642 M35 charge amplifiers should be returned to the manu- 
facturer for calibration. Biased and unbiased responses 
should be checked. Unknown, but apparently unrealistic 
response characteristics made these transducers unusable for 
the reported experiments. 

2 . Piezoresistive Pressure Transducer 

Procurement of a piezoresistive pressure measurement 
system, including a transducer and amplifier such as the 
Kistler Model 603B Series and Model 5004 Charge Amplifier, 
would provided continuous measurements of shock tube 
pressures . 

A second transducer, to be used as a pitot probe for 
flow survey work in shock tunnel applications, is also 
required. Selection would be based on the selected 
nozzle Mach number and probe geometry. Cornell Aeronautical 
Laboratory publications by MacArther should be consulted. 

D. THIN FILM SIGNAL AMPLIFIER 

The current thin film signal amplifier has no record of 
calibration for the past ten years. Voltage fluctuations 
are apparent and dial adjustment response is very poor. 

Prior to the next use the amplifier components should be 
tested and replaced as required. 
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E. DATA ACQUISITION 

It is recommended that the Scientific Alanta signal 
analyzer or its equivalent be used in all future testing. 
Post processing capability, and the availability of four 
channels with independent trigger adjustment provide a 
valuable experimental capability not available with the 
Tektronics 549 and 551 oscilloscopes. 

F. DIAPHRAGMS 

The existing brass diaphragms allow driver operation at 
a constant 600 psi level. Redesign to a higher operational 
level of approximately 700 psi would increase the range of 
reflected pressures available to drive a shock tunnel, with 
an associated cost of higher gas consumption. While not 
essential, this modification would reduce problems which 
would be expected in an extension to a shock tunnel with 
respect to the measurement of small aerodynamic forces. 

Gaydon [Ref. 2:p. 91] and Glass [Ref. 18: pp. 503-511] 
provide valuable design insight into successful diaphragm 
material selection and preparation. 

G. MISCELLANEOUS HARDWARE 

Valve 7 and valve 11 located in the driven end of the 
tube [Ref. 4: Figure 2] exhibited some signs of leakage 
during the loading process. Repair or replacement of the 
valve seats is recommended. 
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APPENDIX B 



ESTIMATION OF SHOCK TUNNEL NOZZLE THROAT AREA 



In order to size a suitable nozzle to be driven by the 
high pressure and temperature in the reflected region, a 
simplified model was developed. The model was based on the 
following assumptions: 

A constant volume reservoir corresponding to the in- 
itial reflected shock conditions was available to drive 
the nozzle flow. The reservoir volume was defined by 
the diameter of the tube and the axial distance from 
the endplate of the tube to the point where the contact 
surface and reflected shock met (Fig. 3); 

The gas remaining in the reflected region expands 
isentropically to fill the constant volume; 

The flow through the nozzle throat was isentropic and 
always choked. (This implies that the evacuated dump 
chamber had sufficient volume, so as not to influence 
the nozzle flow. ) 

The total mass initially in the reflected reservoir is 
given by, 



where 



m. 



R T_ V 0 ' 
g 5 



2 

irD x 



Y 0 = 



(B.l) 



The mass at any later time is given by 



80 



m 



(B . 2 ) 



R T 

g 



t 



V 



0 



As the mass flows through the choked nozzle the rate of 
change of mass in the volume is equal to that which flows 
out through the nozzle, thus 



dm 

dt 




where = 



— 




\R 


2 



V 1 

2(Y 5 -1) 



(B . 3 ) 



Integrating 



4 - f A *P 

- S dm = / _t , dt <B.4> 

o o 

Since the expansion in the reservior is assumed to be isent- 
ropic, the stagnation pressure and temperature are related 
according to 



P 



5 



T 



5 






(B.5) 



Taking the ratio of Equations B.l and B.2 and substituting 
into Equation B.5, 
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1 



m 






Substitution into Equation B.3 gives 



(B . 6 ) 



dm 

dt 




so that 



(B . 7 ) 




Y* + l 



d(-P-) 

m Q 



C 1 A * P 5 
m 0 ^ 0 



t 

/ dt 



( B . 8 ) 



On integrating, the variation of mass in the reflected 
reservior with time is given by 



t 





(B . 9 ) 



and the variation of pressure in the reflected region with 
time is given by 



t 



m 0 /T 5 / 2 \ 

C 1 A * P 5 \ t 5"V 




1 



(B . 10 ) 
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To estimate a suitable nozzle throat, Equation B.10 was 
programmed and a ratio of P t /P 5 equal to one half was 
selected as the criterion for calculating the useful run 
time. Using throat area as a variable all other parameters 
were a function only of the incident shock strength and of 
the physical size of the tube. Thus Equation B.10 enabled 
the available time for steady flow through a given throat 
diameter to be estimated. 
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APPENDIX C 



ESTIMATION OF THE REQUIRED DUMP CHAMBER VOLUME 



In a shock tunnel application, the dump chamber pressure 
is required to remain below the nozzle exit static pressure 
for the useful duration of the ref lected-shock supply condi- 
tions in the shock-tube. To aid in sizing a dump chamber 
for a Naval Postgraduate School shock tunnel the following 
analysis, based on a simplified model, was developed. The 
model assumes the pressure and temperature are uniform in 
the dump chamber, and that the velocity there is effectively 
zero. It is assumed that the flow to the exit of the nozzle 
is isentropic, and that the volume of the nozzle is negligi- 
ble in comparison to that of the chamber. The volume of the 
dump chamber is given in terms of the initial mass contained 
and its properties, by 






RT D m D 



(C.l) 



Dividing by the reflected reservoir volume, similarly writ- 
ten in terms of the initial mass contained and its proper- 
ties, the ratio of dump chamber to reflected shock volume is 
given by 



V 

V 



D 

0 




P D “O 



(C. 2) 
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As the flow continues, the original reflected pressure (P 5 ) 
decreases to lower values (P t ) . Since the nozzle area ratio 
is fixed, the pressure at the exit of the nozzle (P e ) also 
steadily decreases. Meanwhile, the pressure in the dump 
chamber (P D ) steadily increases until, at a point in time, 

P D = P e , and the useful test time is over. The ratio of the 
initial reflected pressure to dump chamber pressure is then 
given by 




P P 
S t 

^p7^ ^p - ^ 

t e 



(C . 3 ) 



where 




A mass balance gives 



m r 



= m Q - 



m . 

i 



+ m 



DO 



(C . 4 ) 



(C . 5 ) 



so that 




1 



m i m D0 
m 0 m 0 



Using Equation B.6 from Appendix B 

J 



(C . 6 ) 
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(C . 7 ) 



m . 

1 



P t 



and substituting in Equation C. 6 , 



1/y 

= i _ / ^t\ 5 + m D0 






P 5 J 



m. 



(C . 8 ) 



Using Equation C.2 applied to the initial conditions, 



m 



DO 



m, 



V 

V, 



D 









(C . 9 ) 



Combining Equations C.2 through C.9 



V, 



D 



V„ 



1 + 



DO 



)©©-< 






D 





M 2 

e a 


Y 5- 1 


1 2 




(*J 


w 




V P 5. 


I 5 





(C.10) 



Equation C.10 provides a relationship between V D /V 0 and the 
value of P t /P 5 which terminates the useful test time, in 
terms of given initial reflected reservoir and dump chamber 
conditions. Using an interval halving technique the dump 
chamber volume requirements were investigated using an 
assumption that the temperature in the dump chamber was 
equal to its initial value, and also with the worst-case 
assumption that it was equal to T 5 . (It is noted that the 
temperature change can be included explicitly through the 
conservation of energy. However the complexity of doing so 
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is not warranted when an estimate of the needed volume is 



all that is required.) 
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APPENDIX D 



SUMMARY OF EXPERIMENTAL DATA 

This appendix contains data taken using thin film gauges 
and the pressure transducer. First, in Figure 42, the 
measurements obtained of shock Mach number using oscillo- 
scopes (preceeding Run 50) are shown. For all runs from Run 
50 on, data were taken using the Scientific Atlanta SD380 
signal analyzer. Tables 6 and 7 list data for runs using a 
mixture of helium and nitrogen and for runs using pure 
helium as a driver, respectively. For each run, trace A 
recorded the response of the thin film at location A, trace 
B recorded the response of the thin film response at station 
B, trace C recorded the response of the endplate pressure 
transducer. Delay times were programmed prior to each run. 
Traces B and C were delayed by the same time. Figure 43 
provides definitions of items which appear on the SD380 
screen and identifies the menu selections of interest. 

The last 35 pages of the appendix contain printouts of 
the SD380 screen for the runs listed in Tables 6 and 7. 

Each printout can be interpreted with reference to Figure 
43. Times were obtained for shock travel between trans- 
ducers on traces A, B and C, using the time scale given on 
the printout, and the delay times listed in Tables 6 and 7. 
These results are included in the two tables. 
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Figure 42. shock Strength vs. P,/P, for Data Reduced 
from Polaroid 
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TABLE 6 



TIME MEASUREMENT 





(70% 


HELIUM 


30% NITROGEN 


DRIVER/N; 


2 DRIVEN) 






Tempera- 


Delay 


Delay 


Time A 


Time B 


Time C 


Run 


ture 
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B 


to B 


to C 


to B 


# 


(F) 


(msec) 


(msec) 


(msec) 


(msec) 


(msec) 


62 


63 


-1.01 


0.00 


2 . 84 


1.67 


3 . 63 


63 
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-1.01 
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1.72 
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-1 . 01 
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70 
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Time C 
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51 
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52 
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SETUP GRP TIME TRC C VW 80DB CH ABCD FR 20KHZ 

01:28:20 SYNC TIME C INP DG X2 WTG H C 2V 
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ETUP GRP TIME QUAD VW 80DB CH ABCD FR 20KHZ 
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SETUP GPP TIME QUAD VW 30DB CH ABCD FP 20KHZ 
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ETUP GRP TIME TRC C VW 80DB CH ABCD FR 20KHZ 

6: 06 SYNC TIME C IMP DG X2 WTG H C 2V 
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GF1P TIME QUAD VW 80DB CH ABCD FR 20KHZ 
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SETUP GRP TIME DUAL VW 80DB CH A0 FR 20KH 
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SETUP GRP TIME DUAL VW 80CB CH AB FR 20KHZ 

00: 33: 37 SYNC TIME A IMP DG X2 WTG H A 5V B 
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APPENDIX E 



COMPUTER PROGRAM FOR SHOCK TUNNEL EVALUATION 

The attached program was used for the evaluation of the 
NPS shock tube and shock tunnel. Sources of shock tube 
relationships are noted in the program were applicable. 
Output was tabulated by the program and graphical analysis 
was provided by the plotting routine ' EASYPLOT' installed on 
the NPS IBM computer. 
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FILE: SH0CK2 FORTRAN A1 



CC THIS PROGRAM CALCULATES THE IDEAL PERFORMANCE OF A SHOCK 

CC TUBE USING THE ANALYSIS PRESENT ED IN THE ’FUNDAMENTALS 

CC OF GAS DYNAMICS’ BY J . A. OHCZAREK PG 395-398 (FIG 8-22) 

8 ’THE SHOCK TUBE IN HIGH TEMPERATURE CHEMICAL PHYSICS’ BY GRA 
GRAYDON HURLE 1963. 



IMPLICIT REAL x 4 (A-H,N,P-Z) 

REAL FCN,X1,X2,XT0L,FTGL, FCN1, FCN2, FCN3, FCN4 
INTEGER I, NLIM,J, 11,112 
EXTERNAL FCN, FCN1 , FCN2 , FCN3 , FCN4 

DIMENSION M 1(200), M2(200), M3 (200), P4P1(200), P2P 1(200) 

DIMENSION P3P4(200) , MS (200) ,US(200), T2 ( 20 0 ) , P5P2 ( 20 0 ) , P5P4(200) 
DIMENSION C3(200),C4(200),P4P3(200),T3(200) , TC ( 20 0 ) , XC ( 20 0 ) 
DIMENSION VC (200) ,P5P1( 200) ,T5T 1(200 ), UR (200), DEL TIM( 200) 
DIMENSION DIST(200),TCC(200),XC1(200), DEL 1 A ( 20 0 ) , PDI ( 8 ) 

DIMENSION XMS(200),XP2P1(200),XP5P2(200),XP3P4(200),XR02R1(200) 
DIMENSION XT2T1(200),URUS(200),TSTAP> 7(200), VDV 0(8), VD(8),AA(8) 
DIMENSION P5(200 ),T5( 200 ),DTHR(200 ),DK 200),MT(200),DTEST( 200) 
DIMENSION PT(200),TT(200),T(200, A ) , TMASS( 200 ) , DT DTHR( 5 ) , PI ( 200 ) 
DIMENSION FT ( A , 1 0 ) , PTPTO ( 1 0 ) 

COMMON /B/ R1 , MM, C8 , DTHR , T 

COMMON /A/ CONST, GAMMA 1, GAMMA4, AAAI SQ, Cl ,C1SQ,CC3, CCA, C5, AAAI ,C9 

COMMON /C/ ATASTR, MT , DT EST 

COMMON /E/ DTHR,T,MT,DTHRDT,MM,C8,R1 

COMMON /D/ P5, T5,T1,AA, PDI, J, 11,112 

COMMON /F/ PTPTO , FT 

DATA XI ,X2, XTOL, FTOL , I , NLIM/0 .1,6.0,0.0001,0.00001,0,50/ 



SELECT GAS CONSTANTS 



GAS 

AIR 

ARGON 

HE 

N2 

C02 



MM 

28.97 
39. 94 
A 

28 . 02 
44 .01 



GAMMA 


SPEED OF 


1.4 


345 M/S 


1 .67 


321 


1.67 


1062 


1 .4 


350 


1 .29 


269 



DRIVEN GAS CONSTANTS 
Rl=8314/28 . 02 
T 1 =288 
GAMMA 1 = 1 .4 
PI =1 . 0E04 

DRIVER GAS CONSTANTS 

R4=831 4/28 .02 

T4 = 288 

GAMMA4 = 1 .40 

P4 = 6 0 0X6 . 8 948 EO 3 

A IS INCREMENT ON P4P1 LOOP 

A= 1 0 0 

MM IS THE NUMBER OF LOOPS 
MM = 6 

P4P1 ( 1 ) =20 



INPUT TUNNEL PARAMETERS( DRIVER LENGTH LI & L4) 
Ll=6 .453 
L 4 = 2 . 37 5 
DT 1 =3* . 0 254 

C TO EVALUATE TEST CONDITIONS SET FL2=1 

THR1 = .0065 
FL 1 = 0 

C FL 2 CALCULATES TEST SECTION MACH # 

FL 2 = 1 

C FL 3 IS FOR NOZZLE SIZING CALCULATIONS 

FL 3 = 0 

C CALCULATES TEST SECTION CONDITIONS 

FL 4 = 1 
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SH0CK2 FORTRAN A1 



FILE: 



C FL6 EXPERIMENT OUTPUT RUNS 

FL6 = 0 

C FL7 = 1 PRINTS MISC RELATIONSHIPS 

FL7 = 0 

C FL8 = 1 EVALUATES PT/PTO RATIOS FOR GIVEN THROAT 

FL8 = 1 

C EVALUATE SPEED OF SOUND CONSTANTS 

A1=SQRT(GAMMA1XR1XT1 ) 

AA=SQRT (GAMMAA*RA*TA) 

AAAI =AA/ A1 
AAA 1 SQ = AA Al**2 



EVALUATE CONSTANTS USED COMMONLY USED 
Cl =( GAMMA 1+1 )/ ( 2XGAMMA 1 ) 

Cl SQ r Cl **2 

C5= C -2XGAMMAA )/ C GAMMA A- 1 ) 

C6 - C GAMMAA-1 )/ C 2XGAMMA A ) 

C7 = C GAMMAA + 1 )/ C AXGAMMAA ) 

C8= C GAMMA1 - 1 )/ C 2XGAMMA 1 ) X ( 1 ) 

C9=( GAMMA 1+1 )/(2X(GAMMAl-l) ) 

Cl 0 = ( GAMMA1/ ( GAMMA1 -1 ) ) 



THIS LOOP INCREMENTS PRESSURE PA/PI BY A AND THEN 
USES INTERVAL HALVING TO DETERMINE M3. 

DO 1 J=1,MM,1 
CONST = P AP 1 ( J ) 

CALL INTHVC FCN,X1,X2,XR,XT0L, FTOL,NLIM, 1 ) 

M3( J)=XR 
C3C J ) =CC3 
CA C J ) = CCA 
XI = XR 
X2 = 6 . 0 
1 = 0 

THE ROOT HAS BEEN FOUND AT XR IF I EQUALS 1 OR 2 
K = J + 1 

PAP1(K)=PAP1( J)+A 
1 CONTINUE 

THIS LOOP USES PA/PI 8 M3 TO DETERMINE P2,P1,MS,M2,U5,P3PA 

DO 5 1=1/ MM , 1 

FIND P2/P1 ( EQN 8-57 ) 

P2P1(I)=PAP1(I)/CA(I) 

M2(I )=SQRT(2X(P2P1(I )+ ( 1/P2P1 ( I ) ) -2 )/ ( GAMMA1X( GAMMA1-1 ) X 
8 C ( GAMMA1 + 1 )/ ( GAMMA1 -1 ) + P2P1 ( I)) ) ) 

FIND MS 8 SPEED OF SHOCK (EQN (8-51) 

MS (I )=SQRT( (P2P1(I)-1)XC1+1) 

USC I ) =MS( I ) XA1 
FIND P3/PA 
P 3PA(I)=C3(I)**C5 
PAP3(I)=1/P3PA( I) 

EVALUATE MAXIMUM DURATION OF UNIFORM FLOW IN THE REGION 
BETWEEN THE SHOCK WAVE AND THE CONTACT SURFACE. 

REFERENCE PAGE A33 PROBLEM 

T2(I)=(2/(P3PA(I)x*C7)-(A*AA*(l-P3PA(I)x*C6)*(l/P3PA(I))x*C7)/ 
&( (GAMMAA-1 )XA1XSQRT( 1+(GAMMA1+1)X(P2P1( I )-l )/ 
S(2XGAMMA1))))XLA/AA 



T 3 IS THE THEORITI CA L MAXIMUM DURATION OF UNIFORM FLOW 
IN THE REGION BETWEEN THE CONTACT SURFACE AND THE TAIL OF 
CENTERED EXPANSION WAVE (REGION 3). ( PG A33 PBLM 8-6) 

T3( I )=( LA/AA)X( (GAMMAA-1 )/2)x( (PAP3( I )xx( 3-GAMMA A) ) 

8XX( 1/ ( GAMMA Ax A) ) )/( l-( P3PA( I )XXC6 ) ) 



TC 8 XC ARE THE TIME AND DISTANCE AT WHICH THE WAVE FORMED BY 
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FILE: SH0CK2 



FORTRAN A 1 



THE REFLECTION OF THE CENTERED EXPANSION WAVE OVERTAKES THE 
CONTACT SURFACE. (REF PG433 PBLM 8-7) 

TC(I)=(L4/A4)X2*(P4P3(I)XXC7) 

XC(I) = (L4X<4/(GAMMAA-1)X(1-(P3P4(I)*XC6)/(P4P3(I)XXC7))) 



CALCULATION OF REFLECTED REGION PARAMETERS ACCORDING TO REF 2 

VELOCITY OF CONTACT SURFACE VC 

VCC I)=2*Al/( GAMMA 1+1 ) X ( MS ( I ) - 1/MS ( I ) ) 

PRESSURE AND TEMPERATURE IN REFLECTED REGION 5 ( EQN 2.31 ; 
P5P1 ( I ) = ( ( 2XGAMMA1 XMS ( I ) XX2- ( GAMMA 1 -1 ) )/(GAMMAl+l ) ) 
8X(( (3XGAMMA 1-1 )XMS( I) XX2-2XC GAMMA 1 -!))/( (GAMMA 1-1) XMS( I )x^2 + 2)) 
T5T 1 ( I) = ( (2X( GAMMA 1-1 )XMS(I )x*^ ( 3-GAMMA 1 ) ) 

8X( (3XGAMMA1-1)XMS( I)XX2-2X(GAMMA1-1 ) )/ ( ( GAMMA 1 + 1 ) XX2XMS ( I ) xx 2 ) ) 
CALCULATE VELOCITY OF REFLECTED SHOCK UR (EQN 2.38) 
UR(I)=US(I)X( (2+(2/(GAMMAl-l))X(l/P2Pl(I)))/ 

&( (GAMMA1 + 1 )/(GAMMAl-l )-l/P2Pl(I )) ) 

OBSERVATION TIME AND DISTANCE I NREFL ECT ED SHOCK 

DELTIM(I)=L1X(US(I)-VC(I))/(US(I)X(UR(I)+VC(I )) ) 
DIST(I)=DELTIM(I)XUR(I) 

TIME TO EXPANSION WAVE REACH CONTACT SURFACE (EQN A. 7) 

C=( (GAMMA A+l )/(2X( GAMMA A-l) ) ) 

TCC( I )=2XLAX(1+(GAMMAA-1 )XM3(I )/2)xxc/AA 
XC1(I)=TCC(I)X( (2XA1X(MS(I)XX2-1))/((GAMMA1+1 )XMS(I))) 
TEST TIME BETWEEN CONTACT SURFACE AND SHOCK 
DELTA(I)= TCC(I)x(l-2X(MS(I)xX2-l)/( (GAMMA1+1 )XMS(I )XX2)) 
P5P2(I)=P5P1(I)/P2P1(I) 

5 CONTINUE 



SECOND METHOD OF CALCULATIONS 
Xl = .l 
X2 = 6 . 0 
XTOL = 0 . 0001 
FTOL=0. 00001 
1 = 0 

NL IM = 5 0 
XR = . 5 

DO 2 J=1 , MM, 1 
CONST = P API ( J ) 

CALL INTHV(FCN1,X1,X2,XR,XT0L,FT0L,NLIM,I) 

XMS ( J ) - X R 

XP2PK J ) = (2XGAMMA1XXMS( J ) *X2-( GAMMA1 -1 ))/( GAMMA 1+1 ) 

XP3PAC J )=1/(PAP1 ( J )/XP2Pl ( J ) ) 

XR02RK J)=( (GAMMA1+1 )XXMS( J)XX2)/( (GAMMA1-1 )XXMS( J )^X2+2) 
XT2T1( J) = ( ( GAMMA 1XXMS( J ) xx2-( GAMMA1 -1 )/2 ) )X( (GAMMA 1-1 )/ 2 XXMS( J)XX 
82+1 )/ ( ( ( GAMMA1+1 )/ 2 ) XX2XXMS( J ) XX2 ) 

XP2P1( J) = (2XGAMMA1XMS( J )XX2- (GAMMA 1-1 ) )/( GAMMA 1 + 1 ) 

XP5P2C J) = ( ( (GAMMA1 + 1)/ (GAMMA 1-1 ) ) + 2-l/XP2Pl ( J ) )/ 

&( l + ( GAMMA1 + 1 )/( GAMMA 1 -1) X1/XP2PK J) ) 

XP3PAC J)=1/(PAP1( J)/XP2P1( J)) 

XR02R1(J)=( (GAMMA 1+1 )XMS(J)XX2)/( (GAMMA1-1 )XMS(J)XX2+2) 
XT2T1(J) = ( ( GAMMA 1 XMS ( J )XX2-( GAMMA 1-1 )/ 2 ) )X( (GAMMA 1-1 )/2*MS(J >xx 
82+1 )/(( (GAMMA 1 + 1 )/ 2 ) XX2XMS ( J ) XX2 ) 

XI =XR 
X2 = 6 . 0 
1 = 0 

M3 HAS BEEN DETERMINED BY INTERVAL HALVING 
DATA XI, X2, XTOL, FTOL,l,N LI M/XR, 5. 0,0. 0001, 0.00001, 0,50/ 

URUS( J)=UR(J)/US( J) 

P5P4( J)=P5P1( J)/PAP1( J) 

C P5( J ) =P5P1 ( J ) XP1 

P5 ( J ) =P5P A ( J ) xpA 
Pl( J)=PA/PAP1( J) 

T5( J)=T5T1( J)XT1 

W R I T E ( X , x ) P5P1( J),P4P1( J),P1( J),UR( J) ,DELTIM( J) 

2 CONTINUE 



C 
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C = =============================================================== 

IF ( FL 3 . EQ . 0 ) GO TO 36 

C THIS SECTION DETERMINES P5 T5 AND DETERMINES TESTIME AS A 

C FUNCTION OF NOZZLE THROAT DIAMETER. 

NN = 1 

DO 3 1=1, MM 

C REFLECTED MO AVAILABLE TO DRIVE NOZZLE 

C=SQRT ( GAMMA 1/R1 )*C 1-C GAMMA1-1 )/2 ) ** ( -1*C9 ) 

TMASSC I)=P5C I )*(URC I )*DELTIMC I ) )*3 . 1A17*DT1**2/CR1*T5(I)*A) 

C WRITEC*,*) P5P1 ( I ) , P5C I ) , T5C I ) , NN, TMASSC I ) 

CALL tTOZZL E( P5( I ) , T5 ( I ) , NN, TMASS ( I ) , C, FL 3 ) 

C601 FORMAT ( , PAP1= , ,D9.3,2X, , P1=’,D9.3,2X, ’TOTAL MASS MO= 1 , D9 . A , 2X, 

C & 1 RUN = 1 , 1 5/ ) 

NN=NN+ 1 

3 CONTINUE 

DO 309 1=1,100 

WRI TEC 8,602) DTHR ( I ) , T ( 1 , 1 ) , T ( 1 , 2 ) , T ( 1 , 3 ) , T ( I , A ) 

602 FORMAT (5(E9.3,2X)) 

C20 1 FORMAT ( 1 P5= 1 , E9 . 2 , 2X, 1 T5= 1 , E9 . 2, 2X, , NN= , ,I1 

C & , 2X, ’TOTAL MASS M0= f ,E9.2) 

309 CONTINUE 
C WRIT EC 8,6 03) 

C6 03 FORMAT ( / , ’THROAT DIAMS2X, 1 RUN1 1 ,6X, 'RUNZ’^X, 'RUNS', 7X, ' RUNA 1 , / ) 
C HRI TECS, 6 01) C PAP1CI),P1CI), TMASSC I), I, 1 = 1, A) 

Z --------------------------------------------------------------- - 

C= = = = = =- = = = = - = = = - = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = -- = = = = = = = = = = = = = = 

36 IF ( FL 8 . EQ . 0 )G0 TO 37 

C THIS SECTION DETERMINES P5 T5 AND DETERMINES TESTIME AS A 

C FUNCTION OF PT/PTO FOR A GIVEN THROAT DIAMETER 

NN = 1 

DO 82 1=1, MM 

C REFLECTED MO AVAILABLE TO DRIVE NOZZLE 

C = SQRTC GAMMA 1/R1 ) x( 1-C GAMMA1-1 )/2 )**C -1*C9 ) 

TMASSC I)=P5CI)*CURC I )*DELTIMC I ) ) *3 . 1 A1 7*DT1**2/ C R1*T5C I )*A ) 
CALL N0ZLE2CP5C I),T5( I ) , NN , TMASSC I ) , C, FL 3 ) 

NN=»N+1 
82 CONTINUE 

DO 89 1=1, MM 

WRI TEC 8,802) PTPTO C I ) , FT ( 1 , 1 ) , FT C 1 , 2 ) , FI C 1 , 3 ) , FTC I , A ) 

802 FORMAT (5(E9.3,2X)) 

89 CONTINUE 

WR ITE(8,803) C PA PI C I ) , I = 1 , MM ) 

8 03 FORMAT C/, AC 1 P A P 1 = 1 , E9 . 3 , 2X ) ) 

Z- ----------------------------------------- - - -- -- -- -- ------------ 

37 IF ( FL 2 . EQ . 0 )GO TO 133 

C DETERMINE TEST SECTION MACH fi AS FUNCTION OF THROAT 

C 

C CALL EXPANC THR1 , DTI ) 

C SUBROUTINE EXPANC THR1 , DTI ) 

C REAL XI, X2, XTOL, FTOL , DT1,THR1 

C INTEGER I , NL I M, J 

C DIMENSION MTC200), DTESTC200), DTDTHRC5) 

CC COMMON /A/ CONST, GAMMA 1 , GAMMA A ,AAA1SQ,C1, C1SQ, CC3, CCA , C5, AA A1 
C COMMON /C/ ATASTR, MT, DTEST , C9 , DTDTHR 

WRI TEC 8,300 ) 

DTEST C 1 ) = DTI 
Xl=.10 
XR = 5 0 

DO 30 J = 1 , 8 
X2=l 00 
XTOL = 0 .0001 
FTOL = 0 . 00001 
1=0 

NL IM=50 

ATASTR = ( DTEST C J)/THR1 )XX2 
C WRIT EC X, X) THR1, DTEST (J), ATASTR, C9 

CALL INTHVC FCN3,X1,X2,XR,XT0L, FTOL,NLIM, I) 

MT C J ) =XR 
XI =XR 
X2 = 1 0 0 
C 

K= J + l 130 
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DTESTCK)=CDTESTCJ)+DT1*.5) 

DTDTHRC J)=DTESTC J)/DT1 
PRINT X, 'HELP 1 

MR I TEC 8, 301 ) MTC J ) , DTDTHR C J ) , DTESTC J ) , ATAST R , THR1 , DTI 
30 0 FORMAT (/, AX, 1 MTEST 1 , IX , 1 DT EST/ DTUBE MX , ’ DTEST ’ , 2X, 1 ATASTR 1 
8, AX, f DTHROAT’ ,3X, ’DTUBE’ ) 

301 FORMAT CACF9.3),2C2X, E9 .3)) 

30 CONTINUE 



133 I FC FL A . EQ . 0 ) GO TO 105 

THIS SECTION EVALUATES TEST CONDTIONS FOR VARIOUS 
TEST CHAMBER AREAS GIVEN A FIXED NOZZLE THROAT DIAMETER 
OUTER LOOP IS FOR PAP1 CHANGE, INNER LOOP EVALUATE CH AREA 
DO 22 J = 1 , MM 
MRITEC8, 400) 

DO 21 1=1,8 
MRITEC*,*) CIO 

C MRITEC*, A 05) MT ( I ) , P5 ( J ) , T5C J ) , PT ( I ) , TT ( I ) , Cl 0 

CA 06 FORMAT C2X, ’MACH’,6X, ’P5 = ’ , 6X, ’ PT ’ , 5X, ’ T5 ’ , 8X , ’ TT ’ , 5X , ’ TSTART ’ ) 

CA 05 FORMAT C 6 C D9 . 2 ) ) 

CALL TCONDCMTCI ),P5CJ),T5CJ),PTCI),TTCI),TSTART(I),THR1 ,C10) 
WRIT EC 8, A 01) MTCI),P5C J ) , PTC I ) , T5C J ) , TTC I ) , MSC J ) , PAP1 C J ) , PI C J ) 
&, TSTART Cl) 

A 0 0 FORMATC/, ’TEST MACH ’, 5X, ’ P5 ’, AX, * TEST P * , AX , ’ T5 ’ , AX , 

& * T EST TEMP’ ,3X, ’MS’, AX, ’PAP1’,6X,’P1’,7X, ’ 1 START ’ ) 

A 0 1 FORMATC9CE8 .3,1X),E6 .1) 

21 CONTINUE 

C DUMP CHAMBER SIZING EVALUATION 

MRITEC8, All) 

DO A61 112=1,1 
DO A60 11=1,1 
C AAC I 1) = . 1*1 1 

AAC 1 1 ) = .8 
PDI C I 1 ) = . 013332 
Xl=-10000 
XR = 50 
X2 =100000 
XTOL=.0001 
FT OL = . 000001 
1 = 0 

NLIM=200 

C MRI T EC 8 , * ) AAC 1 1 ) , PDI C 1 1 ) 

CALL INTHVCFCNA,X1,X2,XR,XT0L,FT0L,NLIM, I ) 

VDVOC II )=XR 

C MRI TEC 8,X) AAC ID, PDICI1), VDVOC ID 

XI = XR 
X2 = 1 0 

VDC I1)=VDV0CI1)XCURC J)*DELTIMC J)*3 . 1 A1 7 * DT 1**2/ A ) 

MRI TEC 8, A 12) P5 C J ) , T5 C J ) , V DVO C 1 1) , V DC 1 1) , AA C 1 1) 

All FORMATC/, 3X, ’ P5 ’ , 8X, ’ T5 ’ , 8X, ’ VDVO ’ , 7X, ’ VD ’ , 8X, ’PTPTO’, 

86X ) 

A 1 2 FORMAT C5CE9.2,1X)) 

A60 CONTINUE 

A61 CONTINUE 

MRITEC 8 , A0 2 ) DEL TI MC J ) , T C 1 , J ) 

A 0 2 FORMATC/, ’TEST TIME OF SHOCK TUBE = ’,D9.A,/ 

8, ’TIME FOR . 5*P5 OUT OF NOZZLE =’,D9.A) 

C 8, ’TIME TO START NOZZLE CEMPERICAL ESTIMAT E) ’ , D9 . A ) 

22 CONTINUE 

C== r ================= = ==== = ====================================== 

105 K=MM 

IF CFL6 . EQ . 0 ) GO TO 699 
MRITEC 8 , 6 95 ) 

MR I TEC 8, 6 96 ) C PAP1 C J ) , MSC J ) , USC J ) , VCC J ) , URUSC J ) , T5C J ) 
8,P5CJ),J=1,K) 

695 FORMATC 2X, ’ XP API ’ , 7X, ’ MS ’ , 9X, ’ VS ’ , 9X , ’ VC ’ , 9X , ’URUS’, 

85X, ’ P5 ’ , 8X, ’ T5 ’ ) 

6 96 FORMAT C6CF9.3,2X),E9 . 3) 

C=========== ======== =========================== =================== 

699 IF CFL1 . EQ. 0)GO TO 108 

MRI TEC 8 ,201) ]_3]_ 
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WRIT EC 8, 20 0 ) ( PAP1 C J ) , MSC J ) , XP2P1 C J ) , XP5P2C J ) , XP3PA C J ) 
8,XR02R1( J),XT2T1( J),XMS(J),URUS(J), J=1,K) 

200 FORMAT (9CF11 . 3 , 2X ) ) 

201 FORMAT (2X, 'XPAP1 1 , 7X, T MS ' , 8X, f XP2Pl 1 , 8X, 'XP5P2',6X, 'XP3PA' 
8, AX, ’ R0U2/R0U1 ' , AX, 'XT2/T1’,6X, ’XMS',6X, 1 UR/US f ) 

WRIT E( 8, 203) C PAP1 ( J ) , MSC J ) , XT2T1 ( J ) , T5T1 ( J ) , P2P1 ( J ) 
&,P5PA(J),DELTIM(J), J=1,K) 

203 FORMAT (7(D9.A,2X)) 

108 HR I TE(8,20A) 

WRIT EC 8, 205) R 1 , T 1 , GAMMA1 , R A , T A , GAMMAA 
20 A FORMAT (3X, 'PAP1',8X, 'MS',8X, I XT2T1 , ,7X, 'T5T1',7X, 'P2P1' 

8,7X, ' P5PA 1 , 5X , ' DELTIME' ) 

205 FORMAT ( f Rl= f ,F9.2,2X, 'T1=’,F9.2,2X, f GAMMAl = T , F9 . 2, 

8/ , 'RA=',F9.2,2X, f TA= f ,F9.2,2X, ' GAMMA A= ' , F9 . 2 ) 



IF CFL7.EQ.0) GOTO 109 
C OUTPUT DATA TO A FILE 

WRITEC8, 103) 

WRIT EC 8,100) CPAP1CI ) , MS ( I ) , M2 C I ) ,M3CI),AAA1, 

8P5P1 Cl ), T5T1CI ),P3PA(I ),P2P1CI ) ,P5P2( I), DELTIMCI ), 1=1, K) 
WRITEC8, 106) 

WRI TEC 8,107) (PAP1(I),URCI),DIST(I),XC1CI),TCC (I), DELTIMCI), 
&DELTACI ),T2CI),I=1,K) 

106 FORMAT (AX, T PAP1 T ,2X, ’VEL R EFL T , AX , 1 DI S REFL f , AX, 

8 1 EXP POS * , AX , f EXP TIME 1 ,AX, 'DEL TIM ', 5X, ' DELTA ', 8X, f T2 f ) 

WRIT E(8,10A) 

107 FORMAT ( 8 C F9 . A , 2X) ) 

WRI TEC 8, 105) (PAP1(I),T2(I),T3CI),TC(I),XC(I),I=1,K) 

100 FORMAT ( 1 1 ( F9 . A , IX)) 

102 FORMAT (AX, f PAPl’,8X, 'T2',8X, , P3PA',6X, 'P2P1 T ) 

103 FORMAT (AX, 'PAP1 ' ,6X, MIS 1 , 

89X, ' M 2 ' , 6 X , ' M3 1 , 8X , 'AAAI »,6X, 'P5Pl f 
8 , 6 X, ' T5T 1 1 , 6X , 

8 IX, 'P3PA ' ,6X, 'P2P1 1 ,6X, l P5P2 l , 6X, 'TESTTIME' ) 

1 0 A FORMAT ( 2X, ' PAP1 ' , 1 2X, ' T2 ' , 8X, ' T3 ' , 1 OX, ' TC ' , 9X, ' XC ’ ) 

105 FORMAT (5CF10 .A,2X)) 

109 X=1 

STOP 
EfID 



REAL FUNCTION FCNCX) 

REAL X 

COMMON/ A/ CONST, GAMMA 1, GAMMA A, AAA 1SQ, Cl, Cl SQ,CC3, CCA, C5, AAAI, C9 
CC3=(l+CGAMMAA-l)*X/2) 

CCA=(CC3xx(2XGAf1MAA) )xy( l/( GAMMAA- 1 ) ) 

FCN = ( 1 + C . 5XAAA1SQ)X( ( GAMfTAl XX/CC3 ) **2 ) 

8XCC1+SQRTC (C1SQ + A/C AAAI SQx( GAMMA lxX/CC3)X)f 2) ) ) )) 

8XCCA-C0NST 

WRI TEC 8, 20 0) FCN, GAMMA1 , GAMMAA , AAA1SQ, CC3, CCA, CONST 
200 FORMAT C 7 ( F8 . 2) ) 

RETURN 

END 

REAL FUNCTION FCN1CY) 

REAL Y 

COMMON/ A/ CONST, GAMMA1 , GAMMAA, AAAI SQ, Cl , Cl SQ,CC3, CCA, C5, AAAI, C9 
FCN1 = ( (2*GAMMAl*YXX2-GAMHAl+l)/(GAMMAl+l))X(l-( (GAMMA A-l)/( AAAI 
8X( GAMMA 1+1 ) )X( Y-l/Y) ) )XXC5-C0NST 
PRINTS, FCN1, CONST, AAA1,C5 
R ETURN 
END 



REAL FUfICT I ON FCN3CZ) 

REAL Z 

COMMON /C/ ATASTR , MT, DTEST 

COMMON /A/ CONST , GAMMA 1 , GAMMAA , AAA1SQ , Cl , C1SQ , CC3, CCA , C5 , AAAI , C9 
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C WRITE(*,x) ATASTR, GAMMA 1,Z,C9 

FCH3=ATASTR-1/ZX(2X(1+(GAMMA1-1 )/2*Z**2) 

&/CCGAMMAl+l)))XXC9 
C WRITECX,x) ATASTR, GAMMA1,Z,FCN3 

RETURN 
END 
C 

REAL FUNCTION FCN4CW) 

REAL W 

C INTEGER J,I2,I1 

DIMENSION P5C200) ,T5C200) ,AAC 8 ) ,PDI( 8 ) ,MTC200 ) 

COMMON/C/ ATASTR 7 MT , DTEST 
COMMON/ D/ P5,T5,T1, AA, PDI, J, II, 112 

FCN4=W-(1+PDI(I1)XT5( J)XW/(P5( J)XT1)-AA( Il)xx.714)xTlx 
8(1 + . 2*100 )xx3 . 5/CAACI1 )*T5C J) ) 

C WRIT EC 8, 131 ) AACI1),PDI( II ),P5( J) , W,T5C J), J, II 

C 131 FORMAT ( ' A = f , E9 . 2 , , PDI= , ,E9.2, 1 P5 = 1 , E9 . 2 , 

C & 1 W = 1 , E9 .2, 1 X , E 9 .2, 1X,F4.1,1X,F4.1) 

RETURN 

END 

C 

C 

SUBROUTINE NOZZL EC X, Y, NN , MTOTAL , C0N3T2 , FL 3 ) 

C THIS ROUTINE DETERMINES TIME FOR NOZZLE FLOW FOR A GIVEN THROAT 

C DIAMETER. 

REAL 5<8 X,Y, MTOTAL, C0NST2, PI 
INTEGER NN, FL 3 

DIMENSION DTHRC200),TC200,4) 

COMMON /A/ C0NST,GAMMA1,GAMMA4, A4A1SQ,C1,C1SQ,CC3,CC4,C5, A4A1,C9 

C COMMON /B/ P5,T5,R1,MM,C8,DTHR,DTHRDT,T 

COMMON /B/ R1,MM,C8,DTHR,T 

PI=3 . 1416 
DTHR(1)=. 0065 
DO 30 1=1,100,1 

T( I, NN)=MT0TALXDSQRT(Y)X2XC1/ .5XXC8-1 )/CC0NST2XPIxDTHRCI)xx2/4*X 
&X( GAMMA 1-1 ) ) 

J = I + 1 

DTHRCJ)=DTHRCI)+.0005 
C WRITEC 8,200) DTHR ( I ) , T C I , NN ) 

C20 0 F0RMAT(2(E9.3,2X)) 

C20 1 FORMAT ( f P5= 1 , E9 . 2 , 2X , 1 T 5= 1 , E9 . 2 , 2X , ' NN= 1 , 1 1 
C & , 2X , 1 TOTAL MASS M0 =, ,E9.2) 

30 CONTINUE 

C WRIT EC 8, 201 ) X, Y, NN, MTOTAL 

RETURN 
END 
C 
C 

SUBROUTINE NOZL E2C X,Y, NN, MTOTAL , C0NST2, FL 3) 

C THIS ROUTINE DETERMINES TIME FOR NOZZLE FLOW FOR A GIVEN THROAT 

C DIAMETER WITH VARAT I ON OF PT/PTO 

REAL 5(8 X,Y, MTOTAL, C0NST2, PI 
INTEGER NN , FL 3 

DIMENSION DTHR(200),TC200,4),PTPTO(10),FT(4,10) 

COMMON /A/ CONST , GAMMA 1 , GAMMA 4 ,A4A1SQ,C1,C1SQ,CC3,CC4,C5,A4A1,C9 
C COMMON /B/ P5 , T5 , R1 , MM, C8, DTHR, DTHR DT,T 

COMMON /B/ R1 , MM, C8 , DTHR, T 
COMMON /F/ PTPTO , FT 
PI =3 . 1 416 
DTHR C 1 ) = . 0 06 5 
PTPTO ( 1 ) = . 2 
DO 80 1=1,4 

FTCI,NN)=MT0TALXDSQRTCY)X2X(1/PTPT0(I)xxC8-1) 

&/CC0NST2xPIxDTHRCl)xx2/«xX 

&XCGAMMA1-D) 

J = 1 + 1 

PTPT0(J)=PTPT0(I)+.2 
80 CONTINUE 

RETURN 
END 
C 

C SUBROUTINE EXPANC THR1 , DTI ) 
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REAL X1,X2,XT0L,FT0L,DT1,THR1 
INTEGER I,NLIM,J 

DIMENSION MTC200) , DTESTC 200) , DTDTHRC 5) 

COMMON /A/ CONST, GAMMA1, GAMMAS, ASA1SQ, Cl, Cl SQ,CC3, CCS, C5, ASA1 ,C9 
COMMON /C/ ATASTR,HT,DTEST,C9,DTDTHR 
DT EST C 1 ) = DT1 
Xl = l .0 
XR = 2 

DO 30 J = 1 , 2 
X2 = 1 0 

XTO L = 0 . 001 
FTO L = 0 . 0001 
C 1 = 0 

NL IM=50 

ATASTR=( DTESTC J)/THR1 )**2 
WRITECX, X) THR1 , DTESTC J ) ,ATASTR 

CALL INTHVCFCN3,X1,X2,XR,XT0L,FT0L,NLIM,I) 

MT(I)=XR 
XI =XR 
K = J + l 

DTEST CK)=CDTESTCJ)+DT1) 

DTDTHRC J)=DTEST(J)/DT1 
WRIT EC 8, 301 ) MTC J) , DTDTHRC J) , DTESTC J) 

30 0 FORMAT C2X, MIT EST 1 ,2X, * DT EST/DTUBE 1 , 2X, 1 DTEST* ) 

WRIT EC 8 , 300 ) 

301 FORMAT C 3 C F9 . 2 ) ) 

30 CONTINUE 
RETURN 
END 



SUBROUTINE TCONDCZl , P5, T5 , PT , TT , TIME, THR1 , CIO ) 

REAL *8 P5,T5, TT, TIME, Z1,THR1,B, CIO 

COMMON /A/ CONST, GAMMA 1, GAMMAS, A A A 1 SQ , Cl , Cl SQ , CC3 , CCS , C5 , ASA1,C9 
COMMON /E/ DTRHC200),TC200), MT , DTHRDT , MM, C8 
COMMON /B/ R1 , MM, C8 , DTHR, T 
WRITECX,S03) 

WRIT EC*, SOS) Z1,P5,PT,T5,TT,C10 
NRITEC*,*) CIO 
A5=SQRTCGAMMA1*R1*T5) 

PT=P5/Cl+( GAMMA 1-1 )/2*Zl **2 ) **C1 0 
TT=T5/C1+CGAMMA1-1)/2*Z1XX2) 

B IS THE ANGLE OUT OF THE THROAT 30 DEGREES 
B=TAN(15*3.1S16/180) 

T I M E = 0 

TIME= .76*THR1/C2*B*A5)*C CDSQRTCZ1) )*C1-1/Z1 ) 

&+l/15*C C DSQRT CZ1 ) ) **3-1 )+l/25*C C DSQRT CZ1 ) ) **5-1 ) ) 

A 03 FORMAT C2X, 'MACH 1 ,6 X, f P5 = 1 , 6 X , * PT * , 5X , 1 T5 ' , 8X , * TT 1 , 5X , * T START ' ) 
WRITEC*,SOS) Z1,P5,PT,T5,TT, CIO, TIME 
SOS FORMAT C 7 C D9 . 2 ) ) 

RETURN 

END 



SUBROUTINE INTHVC FCN, XI , X2, XR, XTOL , FTOL , NLIM, I ) 



THIS ROUTINE MAS DEVELOPED BY PROF. MILLER DEPT AERONAUTICS 
SUBROUTINE INTHV : 

THIS SUBROUTINE FINDS THE ROOTS OF A FUNCTION , 
FCX) = 0, BY INTERVAL HALVING . 



PARAMETERS ARE : 

FCN -FUNCTION THAT COMPUTES VALUES FOR FCX). MUST BE DECLARED 

EXTERNAL IN CALLING PROGRAM. IT HAS ONE ARGUMENT, X. 
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XI, X2 -INITIALS VALUES OF X.F(X) MUST CHANGE SIGN AT THESE PTS 

XR -RETURNS THE ROOT TO THE MAIN PROGRAM. 

XT OL , FTOL -TOLERANCE VALUES FOR X,F(X) TO TERMINATE ITERATIONS. 
NLIM -LIMIT TO NUMBER OF ITERATIONS. 

I -A SIGNAL FOR HON ROUTINE TERMINATED. 

1=1 MEETS TOLERANCE FOR X VALUES. 

1=2 MEETS TOLERANCE FORF(X) 

I=-l NLIM EXCEEDED. 

I =-2 F ( XI ) NOT OPPOSITE IN SIGN TO F(X2) 

WHEN THE SUBROUTINE IS CALLED, THE VALUE OF I INDICATES WHETHER TO 
PRINT EACH VALUE OR NOT. 1=0 MEANS PRINT THEM, I. HE. 0 MEANS DON’T. 



REAL FCN,Xl,X2,XR,XTOL,FTOL 
INTEGER NLIM, I, J 
REAL FI , F2 , FR , XERR 



CHECK THAT F(X1) & F(X2) DIFFER IN SIGN. 

FI = FCN(Xl) 

F2= FCNCX2) 

IF (F1XF2 .GT. 0.0) THEN 
I = -2 

PRINT 201 
RETURN 
END IF 



COMPUTE SEQUENCE OF POINTS CONVERGING TO THE ROOT 

DO 20 J = 1 , NL IM 
XR = ( XI +X2 )/ 2 . 0 
FR = FCN(XR) 

XERR = ABS(Xl-X2)/2.0 
IF (I .EQ. 0.0) THEN 
PRINT 199, J , XR , FR 
END IF 



CHECK ON STOPPING CRITERIA 

IF C XERR .LE. XTOL) THEN 
1 = 1 

PRINT 202, J , XR , FR 
RETURN 
END IF 



IF ( ABS(FR) .LE. FTOL ) THEN 
1 = 2 

PRINT 203, J , XR, FR 
RETURN 
END IF 



COMPUTE NEXT X 

IF ( FR^Fl .GT. 0.0 ) THEN 
XI =XR 
FI = FR 
ELSE 
X2 = XR 
F2 = FR 
END IF 

20 CONTINUE 



C 
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WHEN LOOP IS NORMALLY COMPLETED, NL IM IS EXCEEDED 
I=-l 

PRINT 200 , NLIM, XR, FR 
RETURN 



199 FORMAT ( f AT ITERATION 1 , 13, 3X, 1 X= 1 , El 2 . 5 ,4X, 1 F(X) =',E12.5) 

200 FORMAT ( / 1 TOLERANCE NOT MET AFTER 1 , 14, 1 ITERATIONS X = ', 

1 E12.5, 1 F( X ) = 1 , El 2 . 5 ) 

201 FORMAT ( / 1 FUNCTION HAS THE SAME SIGN AT X1&X2 1 ) 

202 FORMAT ( / 1 TOLERANCE MET IN M4,' ITERATIONS X =',E12.5, 

1 1 FCX) = 1 , El 2 . 5 ) 

203 FORMAT ( / 1 F TOLERANCE MET IN 1 , 14 , 9 ITERATIONS X = ',E12.5, 

1 1 FCX) =',E12.5) 

END 
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